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This  paper  deals  with  the  technology  of  burning  natural  gas  with  combustion  air  preheated  to  1300  °C. 
The  objective  of  this  work  is  to  assess  the  potential  of  several  combustion  models  and  NOx  postprocessors 
for  their  abilities  to  predict  NOx  emissions.  The  steady-state  computational  results  have  been  compared 
both  with  in-furnace  measurements  and  with  the  measured  furnace  exit  parameters.  The  following  main 
conclusions  have  been  drawn:  (1)  with  the  exception  of  the  small  region  located  within  the  natural  gas  jet, 
the  computations  resulted  in  good  quality  predictions;  (2)  the  NOx  has  been  formed  in  a  thin  elongated 
region  (flamelet)  located  between  the  natural  gas  jet  and  the  air  jet;  (3)  the  NO  has  been  formed  mainly 
by  the  thermal  path,  and  the  NO-reburning  mechanism  was  of  little  importance;  and  (4)  there  are  strong 
indications  that  the  non-stationary  behavior  of  the  weak  (fuel)  jet  must  be  accounted  for  if  the  predictions 
within  the  fuel  jet  were  to  be  perfected. 


Introduction 

Combustion  of  fuels  with  high-temperature  air 
and  large  quantities  of  flue  gas  is  perhaps  the  most 
rapidly  developing  combustion  technology  of  the 
past  decade  [1-3].  The  technology,  when  applied  to¬ 
gether  with  an  efficient  heat  recovery  system,  offers 
substantial  energy  savings  and  drastic  reductions  in 
C02,  CO,  and  NOv  emissions.  Initially,  the  technol¬ 
ogy  was  named  in  Japan  as  excess  enthalpy  combus¬ 
tion,  the  term  used  by  Hardesty  and  Weinberg  [4] 
to  describe  burners  with  heat  recycling.  Later  on  it 
was  renamed  high-temperature  air  combustion 
(HTAC).  In  Italy,  the  process  has  been  named  mild 
combustion  [5],  while  in  the  United  States,  the  term 
low-NOx  injection  is  used.  In  Germany,  the  process 
has  been  called  flameless  oxidation  [3]. 

There  are  different  methods  of  realization  of  this 
combustion  technology.  However,  two  of  them  seem 
to  be  widely  applied.  The  flameless  oxidation  burner 
[3]  possesses  a  central  jet  of  natural  gas  and  a  num¬ 
ber  (typically  6  or  12)  of  air  jets  supplying  preheated 
combustion  air.  The  combustion  air  jets  are  spaced 
around  the  central  natural  gas  jet.  Before  mixing 
with  the  natural  gas  jet,  the  air  jets  entrain  large 
quantities  of  the  recirculated  combustion  products. 
Thus,  the  fuel  is  oxidized  in  combustion  products 


containing  some  oxygen  (typically  not  more  than 
5%).  In  Japanese  designs  [6],  in  IFRF  designs  [2], 
the  preheated  air  is  supplied  through  a  central  jet, 
whereas  natural  gas  is  provided  by  several  injectors 
located  on  the  circumference  of  the  burner  (see  also 
Ref.  [7]).  The  natural  gas  injectors  are  positioned 
away  from  the  air  jet  so  as  to  inject  the  fuel  into  hot 
combustion  products  that  contain  l%-2%  oxygen. 
The  natural  gas  jets  do  not  mix  with  the  air  stream 
until  further  downstream  in  the  furnace,  and  by  that 
time  (distance),  the  fuel  jets  are  substantially  diluted 
with  combustion  products.  In  both  flameless  oxida¬ 
tion  burners  and  combustion  systems  with  fuel  in¬ 
jectors  positioned  away  from  the  air  jet,  the  fuel  gas 
is  oxidized  in  an  environment  that  contains  a  sub¬ 
stantial  amount  of  inert  gases  (N2,C02,H20)  and 
some,  typically  not  more  than  5%,  oxygen.  In  this 
respect,  both  modes  of  combustion  are  similar.  In 
most  industrial  applications,  burners  are  used  in 
pairs  with  a  firing  cycle  of  20-50  s  due  to  the  built- 
in  heat  regenerators  or  recuperators. 

There  has  been  relatively  little  support  of  the  fun¬ 
damental  combustion  community  to  provide  the  ad¬ 
equate  scientific  insight,  as  pointed  out  in  Refs.  [1,2]. 
Mathematical  modeling  of  flameless  oxidation  has 
received  little  attention.  The  first  modeling  efforts 
originated  from  Japanese  industry  [8,9].  They  show 
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Precombustor 


Fig.  1.  Experimental  furnace. 


that  numerical  codes  based  on  the  mixture-fraction/ 
probability  density  function  (pdf)  approach  are  able 
to  describe  the  main  characteristics  in  terms  of  the 
flow  and  temperature  fields  of  high  air  temperature 
combustion.  However,  they  found  that  the  present 
NOx  models  might  require  improvement  to  describe 
properly  the  NOx  formation  under  low-temperature 
conditions.  In  the  above-cited  work,  the  assessment 
of  the  model  performance  was  made  exclusively  on 
the  basis  of  comparing  the  predictions  with  the  mea¬ 
sured  furnace  exit  values  since  no  in-furnace  data 
were  available.  Furthermore,  although  the  burners 
operated  in  cycles,  the  predictions  were  carried  out 
under  steady-state  conditions. 

More  recently,  Coelho  and  Peters  [10]  carried  out 
numerical  simulations  of  a  mild  combustion  (FLOX) 
burner  that  operated  at  10  kW  thermal  input  with  a 
relatively  low  level  of  air  preheat  (500  °C).  The  ex¬ 
perimental  furnace  was  operated  at  a  rather  low 
temperature  of  1000  °C,  resulting  in  NOx  emissions 
around  10  ppm.  Coelho  and  Peters  [10]  argued  that 
the  steady  flamelet  library  was  unable  to  “correctly 
describe  the  formation  of  NO  since  this  is  a  chemi¬ 
cally  slow  process,  which  is  sensitive  to  transient  ef¬ 
fects”  while  the  unsteady  flamelet  model  was  able  to 
predict  “the  correct  order  of  magnitude  of  NO  emis¬ 
sions.”  Also  in  this  case,  the  conclusions  were  drawn 
on  the  basis  of  the  comparison  of  the  predictions 
with  the  furnace  exit  values  since  limited  in-furnace 
measurements  were  available.  Another  interesting 
paper  was  published  at  the  time  of  preparing  of  this 
paper.  Pasenti  et  al.  [11]  undertook  a  task  of  simu¬ 
lating  a  200  kW  FLOX  burner  operated  under  cyclic 
(unsteady  state)  conditions.  The  furnace  exit  tem¬ 
perature  was  varied  in  the  range  1090-1330  °C  (us¬ 
ing  a  heat  sink)  resulting  in  NOx  emissions  that  in¬ 
creased  (with  the  furnace  exit  temperature)  from  4 
to  40  ppm  (at  3%  02).  The  steady-state  numerical 
simulations  [11]  resulted  in  good  predictions  of  the 
total  radiative  heat  flux.  The  NOx  emissions  were 
substantially  underpredicted  although  the  depen¬ 
dence  (trend)  of  the  NOx  emissions  on  the  furnace 


exit  temperature  was  correctly  represented  in  the 
predictions. 

It  is  apparent  that  there  is  a  lack  of  in-furnace  data 
corresponding  to  flameless  oxidation,  and  conse¬ 
quently  the  mathematical  modeling  is  speculative.  It 
is  extremely  difficult  to  generate  such  data  under 
non-stationary  conditions  that  are  typical  of  indus¬ 
trial  applications  of  flameless  oxidation  combustion. 
To  the  best  of  the  authors’  knowledge,  there  is  only 
one  set  of  in-furnace  data  collected  under  steady- 
state  conditions  that  contains  the  in-furnace  mea¬ 
sured  temperatures,  velocity  field,  chemistry  field 
(02,C0,C02,  H2,  CH4,  NOx,  NO)  and  radiative 
fluxes.  In  1997,  the  IFRF  [2]  generated  such  data. 
We  have  used  [12]  the  IFRF  data  to  develop  and 
validate  models  for  predicting  the  flow  field  (includ¬ 
ing  jet  entrainments),  temperature  field,  and  chem¬ 
istry  field  (02,  C02,  CO,  H20,  unburned  hydrocar¬ 
bons).  In  this  publication,  we  focus  exclusively  on 
predicting  NOx  emissions.  Later  on,  we  will  dem¬ 
onstrate  that  our  steady-state  calculations  result  in 
excellent  NO  predictions  for  conditions  of  high  air 
preheat  (1300  °C)  and  high  furnace  exit  temperature 
(1222  °C)  and  for  a  steady-state  operation  of  a  nat¬ 
ural  gas  burner  at  0.58  MW  fuel  input. 


Objectives 

The  objective  of  our  work  is  to  develop  and  vali¬ 
date  a  mathematical  model  for  predicting  combus¬ 
tion  characteristics,  including  NOx  emissions,  of 
flameless  combustion  of  natural  gas,  since  the  con¬ 
ditions  of  uniform  and  low- temperature  and  low-02 
concentration  are  substantially  different  from  those 
encountered  in  conventional  burner  technology.  We 
wish  to  identify  the  location  of  regions  of  NOx  for¬ 
mation  and  establish  how  important  are  the  NOx  re¬ 
burning  mechanism  and  prompt-NO  formation 
mechanism. 


Experimental 

The  semi-industrial  scale  experiments  [2]  were 
carried  out  in  a  refractory  lined  furnace  with  a  2  X 
2  m  cross-section  and  a  length  of  6.25  m,  (Fig.  1). 

The  furnace  was  equipped  with  only  one  burner 
operated  under  steady-state  conditions.  The  high- 
temperature  regenerator  was  replaced  by  a  precom¬ 
bustor  where  natural  gas  was  burned  with  air  under 
lean  conditions.  Oxygen  was  added  to  the  hot  gases 
leaving  the  precombustor  to  maintain  its  concentra¬ 
tion  at  21%.  Such  a  comburent  simulated  the  air 
preheat  of  the  excess  enthalpy  combustion  process. 
The  experimental  burner  (see  Fig.  2  of  Ref.  [12]) 
was  a  1:1  copy  of  NFK  burner  that  operated  at 
0.58  MW  fuel  input  and  the  combustion  air  was  pre¬ 
heated  to  1300  °C. 
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TABLE  1 

Experimental  conditions 


Flow  rate 
(kg/h) 

Temp. 

(°C) 

Enthalpy 

(MW) 

Composition  %  (vol) 

Natural  gas 

47 

25 

0,58 

CH4  87.8%,  C,H6  4.6%,  C3H8  1.6%,  C4H01  0.5%, 
5.5%  N2;  LCV  =  44.76  MJ/kg 

Comburent 

830 

1300 

0.35 

19.5%  wet  02,  59.1%  wet  N2,  15%  H20, 

6.4%  wet,  C02,  NO  110  ppm  vd 

Furnace  exit  gases 

877 

1220 

0.380 

1.6%  wet  02,  140  ppm  vd  NOx 

The  comburent  was  supplied  through  the  central 
channel  with  an  injection  velocity  of  —85  m/s.  The 
natural  gas  injectors  were  located  280  mm  away 
form  the  burner  centerline.  The  fuel  was  supplied 
through  two  injectors  with  an  injection  velocity  of 
—100  m/s.  The  input-output  conditions  of  the  ex¬ 
periment  are  listed  in  Table  1 . 

During  the  test,  the  flue  gas  temperature  and  flue 
gas  compositions  were  continuously  monitored  using 
a  suction  pyrometer  and  a  gas  sampling  probe.  Un¬ 
der  steady  operating  conditions,  a  complete  set  of 
in-flame  measurements  was  acquired,  including  the 
flow  and  mixing  patterns,  the  total  radiative  heat  flux, 
and  the  combustion  gas  temperatures  and  gas  com¬ 
position.  Detailed  mapping  of  the  turbulent  veloci¬ 
ties  was  carried  out  using  laser  doppler  velocimetry 
(LDV).  The  in-furnace  temperatures  were  measured 
using  a  suction  pyrometer,  while  the  furnace  wall 
temperatures  were  obtained  using  type  B  thermo¬ 
couples.  A  cranked  probe,  operated  under  quench¬ 
ing  rate  of  10'-108  K/s,  was  used  to  sample  the  gas 
which  was  analyzed  for  oxygen,  carbon  monoxide, 
carbon  dioxide,  methane,  hydrogen,  and  nitric  ox¬ 
ides.  An  ellipsoidal  radiometer  was  used  to  measure 
hemispherical  fluxes  at  the  furnace  wall.  The  in-fur- 
nace  measurements,  reported  in  this  paper,  were  ob¬ 
tained  by  traversing  the  furnace  in  the  horizontal 
plane  of  the  two  fuel  injectors  [2]. 

Mathematical  Models 

Combustion  Models 

In  the  present  study,  three  different  combustion 
models  were  applied  to  simulate  natural  gas  com¬ 
bustion  with  high-temperature  air: 

•  An  eddy  break-up  (EBU)-mixed-in  burned  model 

•  An  eddy  dissipation  (EDC)-chemical  equilibrium 
model 

•  A  mixture  fraction/pdf  model  (PDF)  with  equilib¬ 
rium  look-up  tables. 

Due  to  the  lack  of  space,  only  a  short  description 
of  the  applied  models  is  given  here  and  the  reader 


is  requested  to  consult  the  original  papers.  The 
IFRF  EBU  combustion  model  [13]  considers  six 
major  chemical  species:  hydrocarbons  (CxHy),  oxy¬ 
gen,  carbon  monoxide,  carbon  dioxide,  water  vapor, 
and  nitrogen.  The  combustion  is  modeled  with  a 
two-step  reaction  scheme.  The  hydrocarbon  CxHy, 
is  first  oxidized  to  CO  and  H20,  while  intermediate 
CO  oxidizes  to  C02.  The  various  source  and  sink 
terms  entering  the  chemical  species  mass  balance 
equation  are  described  using  the  eddy  break-up 
model  of  Magnussen  and  Hjertager  [14].  The  eddy 
break-up  model  makes  use  of  a  proportionality  con¬ 
stant,  the  mixing  rate  coefficient,  that  is  calculated 
in  each  numerical  cell  as  a  function  of  the  local  Rey¬ 
nolds  number  [13]. 

The  EDC  model  of  Magnussen  [15]  is  derived  us¬ 
ing  a  detailed  description  of  the  dissipative  process 
in  the  flow.  The  model  derivation  starts  with  the  ex¬ 
pressions  describing  the  dissipation  of  turbulent  ki¬ 
netic  energy.  Further  on,  characteristic  velocity  and 
length  scales  for  turbulent  structures  where  the  mix¬ 
ing  is  completed  are  derived.  In  the  EDC,  these  tur¬ 
bulent  structures  are  called  fine  structures.  Using 
the  velocity  scale  and  length  scale  associated  with 
the  fine  structures,  an  expression  for  calculating  the 
mass  transfer  rate  between  the  surrounding  and  the 
fine  structure  is  derived.  The  chemical  composition 
within  the  fine  structures  is  obtained  by  solving  a 
perfectly  stirred  reactor  problem  using  the  thermo¬ 
dynamic  equilibrium  approach.  The  equilibrium 
composition  was  calculated  using  the  CEA  code  of 
NASA  and  13  species  were  considered:  C2H6,  CH4, 
C02,  02,  H20,  H2,  H,  OH,  O,  NO,  N02,  and  N2. 
The  model  used  here  is  identical  to  the  one  used  by 
Brink  et  al.  [16]  and  Orsino  at  al.  [12]. 

The  third  model  is  the  mixture  fraction  model 
combined  with  the  /^-probability  density  function  to 
handle  the  fluctuations.  The  /^-function  PDF  shape 
is  calculated  from  the  mean  mixture  fraction  and  its 
variance.  An  additional  transport  equation  for  the 
variance  of  the  mixture  fraction  is  solved.  The  com¬ 
putational  time  is  optimized  using  a  preprocessor  to 
obtain  the  equilibrium  composition  of  all  the  chem¬ 
ical  species  considered  in  the  calculations.  During 
the  preprocessing  three-dimensional  look-up  tables 
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are  generated  for  the  mean  flame  temperature,  den¬ 
sity,  and  species  mole  fraction  as  a  function  of  the 
mean  mixture  fraction,  its  variance,  and  the  enthalpy. 
In  the  present  work,  11  chemical  species  (CH4, 

^2^6?  C3H85  to?  co,  co2,  h2o,  oh,  n2,  o2, 

and  H2)  have  been  considered.  While  for  the  EBU 
and  the  EDC,  the  IFRF  software  was  used,  the 
Fluent  CFD  code  (version  5.4)  was  used  to  run  the 
PDF  model. 


NOx  Postprocessors 

Chemical  reaction  rates  of  NO  formation  and  NO 
concentrations  are  postcalculated  using  previously 
computed  velocities,  turbulence,  temperature,  and 
chemistry  fields.  In  the  combustion  of  natural  gases, 
the  two  principal  formation  mechanisms  of  NO  are 
oxidation  of  molecular  nitrogen  (thermal-NO  for¬ 
mation)  and  the  mechanism  in  which  hydrocarbon 
radicals  attack  molecular  nitrogen  to  form  nitrogen 
radicals,  the  latter  subsequently  reacting  with  the 
molecular  oxygen  to  form  NO  (prompt-NO  forma¬ 
tion).  In  fuel-rich  regions,  NO  can  be  reduced  by 
CH,  radicals  to  form  species  like  HCN,  which  sub¬ 
sequently  produces  NO  or  N2.  This  reduction  mech¬ 
anism  is  known  as  NO-reburning  mechanism  or 
fuel-staging  mechanism. 

IFRF  NOx  model 

Thermal  NO  and  prompt  NO  are  calculated  in  the 
same  way  as  in  Ref.  [13]  with  no  changes  to  any  of 
the  model  parameters  and  constants.  However,  we 
added  the  NO-reburning  path  following  the  work  of 
Chen  [17]  and  De  Soete  [18]: 

NO  +  CrHv  ->  HCN  +  .  .  .  (R9) 

HCN  +  02  ->  NO  +  .  .  .  (RIO) 

HCN  +  NO  -►  N2  +  .  .  .  (Rll) 

It  might  be  useful  to  note  that  in  light  of  the  present 

knowledge  of  NO  precursors  kinetics,  the  famous 
CH  +  N2  reaction  does  not  go  to  NCN  and  N  as 
products,  but  rather  the  NCN  +  H  reaction.  Hence, 
the  sequence  of  reactions  R9  +  RIO,  proposed  by 
Chen  [17],  and  used  in  this  study,  does  not  have  a 
firm  basis  in  kinetics,  as  we  now  understand  them. 

For  reaction  R9,  Chen  [17]  proposes  the  rate 

r9_NO  =  2.68  *  106  [Cx  Hy]  [NO] 

18,900\ 

RT  ) 

For  reactions  10  and  11,  De  Soete  [18]  proposes 


r10-NO  =  1010  (—V  [HCN]  [02]6 

\  P  mix  ■ 

(  67,000\ 

exp  r-Rr-j 

ru_NO  =  3  •  1012  [HCN]  [NO] 

'  Pmix  > 

(  60,000\ 
exp  r-Rr-j 

where  Mmix  and  pmix  are  the  mixture  molecular 
weight  and  the  mixture  density,  respectively.  The 
power  b  is  locally  calculated  according  to  De  Soete 
[18]: 

b  =  -1.4799  -  1.4593c  -  0.4859c2 

-  0.0546c3,  c  =  In  [02]j 

The  interaction  between  turbulent  flow  and  chem¬ 
istry  is  taken  into  account  through  the  transport 
equation  of  NO  chemical  species.  The  source  term 
of  this  equation,  which  is  the  time-mean  chemical 
reaction  rates  of  NO  formation/reduction,  is  com¬ 
puted  on  the  basis  of  the  single-variable  pdf  model. 
Since  the  formation  and  reduction  mechanisms  of 
NO  are  strongly  temperature-dependent,  a  pre¬ 
sumed,  single-variable  pdf  approach  has  been  cho¬ 
sen.  For  the  instantaneous  gaseous  phase  tempera¬ 
ture,  T,  the  ft- pdf  /Ljf  (T;  a,b)  has  been  chosen.  The 
time-mean  rates  of  thermal  NO,  prompt  NO,  and 
NO  reburning  are  approximated  by  the  expectations 
Exp(r/,  _  NO)  of  the  original  chemical  reaction  rates 

rk  -  NO 

rTb 

rk  -  no  =  EXP  (rj,  _  NO)  =  rk  -  no  (E)/?pdf 

jk u 

dT 

\T;  a,b)  — - — ,  k  =  t,  p,  r 

-U  )  ~~  h 

Nitric  oxide  models  used  with  the  pdf  combustion 
model 

The  NOx  postprocessing  used  with  the  pdf/com- 
bustion  models  considers  the  thermal  and  prompt 
mechanisms.  The  NOx-reburning  mechanism  is  ex¬ 
cluded.  The  extended  Zeldovich  mechanism  is  used, 
and  the  rate  constants  are  taken  from  Hanson  and 
Salimian  [19]  The  prompt  NOx  mechanism  is  mod¬ 
eled  following  the  overall  simplified  NO-formation 
rate  proposed  by  the  De  Soete  [18].  The  interactions 
with  the  turbulence  are  predicted  using  a  pdf  ap¬ 
proach  using  a  /^-function.  The  limits  of  the  pdf  in¬ 
tegration  are  determined  from  the  maximum  and 
minimum  values  of  the  predicted  temperature  in  the 
computational  domain.  The  Fluent  CFD  code  (ver¬ 
sion  5.4)  is  used  in  this  paper  to  run  the  NOx  post¬ 
processor  described  in  this  paragraph. 
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x  =  0.15  m  x  -  0.73  m  x  =  2.05  m 


Methane  (%,  dry) 


■  IFRF-EBU  °  Fluent5-pdfll 

IFRF-EDC  A  measured 

♦  IFRF-EDC-more-turb 


Fig.  2.  Measured  and  predicted  methane  concentrations. 


Results 

The  experimental  furnace  had  two  symmetry 
planes;  therefore,  a  grid  reproducing  only  one-quar¬ 
ter  of  the  furnace  was  built.  The  grid  was  a  struc¬ 
tured  hexahedral  mesh  of  the  geometrical  domain 
with  92  X  42  X  129  (498,  456)  cells.  The  flow  field 
has  been  calculated  using  the  k-e  model  of  turbu¬ 
lence.  The  furnace  exit  parameters  have  been  pre¬ 
dicted  with  a  very  good  accuracy.  For  example,  a 
temperature  of  1493  K  was  measured  while  the  pre¬ 
dicted  values  were  1522,  1531,  and  1521  K  for  the 
EBU,  Fluent5-pdf,  and  EDC  models,  respectively. 
A  concentration  of  2.5%  02  was  calculated  for  all  the 
models. 

Prior  to  carrying  out  an  assessment  of  the  NOx 
processors,  it  is  imperative  to  identify  the  imperfec¬ 
tions  in  the  computed  in-furnace  temperature  and 
oxygen.  Fig.  2  shows  the  unburned  methane  con¬ 
centrations,  while  the  measured  and  predicted  time 
mean  temperatures  are  shown  in  Fig.  3.  In  these 
figures,  zero  radial  distance  represents  the  furnace 
centerline  (comburent  injection),  and  the  natural  gas 
injection  is  located  at  0.28  m  radial  distance.  For  the 
time  being,  we  shall  concentrate  on  the  predictions 
obtained  with  the  EBU  model  (marked  in  Figs.  2 
and  3  as  IFRF-EBU),  the  EDC  model  (marked  as 
IFRF-EDC),  and  the  pdf  model  (marked  as  Flu- 
ent5-pdfll).  Later  on,  we  will  consider  the  predic¬ 
tions  marked  in  Figs.  2  and  3  as  IFRF-EDC-more- 
turb.  The  figures  demonstrate  the  following 
(predictions  of  the  velocities,  02,  CO,  radiative  heat 


fluxes,  and  the  furnace  wall  temperature  can  be 
found  in  Ref.  [12]): 

1.  The  numerical  models  have  correctly  reproduced 
the  uniformity  of  the  temperature  field  and  the 
oxygen  field.  All  the  three  combustion  models 
provided  the  predictions  of  similar  quality.  The 
EDC  model  and  the  pdf/ mixture  fraction  model 
when  coupled  with  chemical  equilibrium  have 
provided  very  similar  results. 

2.  The  three  tested  combustion  models  cannot  de¬ 
scribe  the  chemistry  and  temperature  field  in  the 
fuel  jet  region. 

The  temperatures  in  the  fuel  jet  are  determined 
by  the  rate  of  energy  release  (due  to  combustion) 
and  the  rate  of  entrainment  of  hot  combustion  prod¬ 
ucts  into  the  jet.  The  oxygen  needed  for  methane 
combustion  is  entrained  into  the  jet  from  the  sur¬ 
roundings.  All  three  models  underestimated  the  rate 
of  methane  consumption  as  shown  in  Fig.  2.  The 
predicted  entrainment  into  the  fuel  jet  was  by  a  fac¬ 
tor  of  2  lower  than  the  measured  one  (computational 
runs  with  the  realizable  k-e  turbulence  model  or  the 
Reynolds  stress  model  offered  marginal  changes). 
Our  analysis  of  jet  entrainment  [12]  has  shown  that 
even  correcting  for  the  entrainment  would  not  result 
in  the  temperature  increase  that  would  fully  com¬ 
pensate  for  such  a  large  difference  between  the  mea¬ 
sured  and  the  predicted  temperatures  within  the 
methane  jet.  The  questions  we  wish  to  address  are: 
What  is  the  reason  for  the  failure  of  all  the  models 
in  the  fuel  jet  region?  What  are  the  consequences  of 
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x  =  0.15  m  x  =  0.43  m  x  =  0.73  m  x  =  1.03  m 


Temperature  (K) 


■  IFRF-EBU  °  Fluent5-pdfll 
IFRF-EDC  A  measured 
*  IFRF-EDC-more-turb 


Fig.  3.  Measured  and  predicted  temperatures. 


TABLE  2 

Measured  and  predicted  NOx  emissions  at  the  furnace  exit 


Combustion  model 

NOx  model 

TotaP 

Predictions  (ppm,  dry) 

Thermal1 

Prompta 

IFRF-EBU 

IFRF 

137 

130 

9 

IFRF-EDC 

IFRF 

132 

130 

6 

Fluent5-pdf  (temp) 

Fluent5 

150 

147 

7 

Fluent5-pdf  (mixture  fraction) 

Fluent5 

130 

127 

5 

IFRF-EDC-more  turbulence 

IFRF 

132 

130 

6 

Note:  Measurements:  Inlet  NOx  (in  comburent),  110  ppm,  dry;  furnace  exit  NOx,  140  ppm,  dry. 
a  In  column  Total,  the  predicted  furnace  NO  concentration  is  listed.  The  predictions  are  obtained  with  both  thermal- 
NO  and  prompt-NO  submodels  activated  and  the  inlet  NO  is  110  ppm,  dry.  Column  Thermal  lists  the  predictions  when 
only  the  thermal-NO  submodel  is  activated  and  the  inlet  NO  (in  comburent)  is  110  ppm,  dry.  Column  Prompt  refers  to 
the  computational  run  with  only  prompt-NO  submodel  activated  and  nil  NO  inlet  concentration. 


the  failure  with  respect  to  the  NOx  emissions  pre¬ 
dictions?  We  begin  with  the  second  question. 

NOx  Etnissions 

The  same  NOx  model  (described  above)  was  used 
with  the  EDC  and  EBU  combustion  model.  Thus, 
the  differences  in  the  NOx  predictions  between 
these  two  cases  are  solely  due  to  differences  in  the 
temperature  and  chemistry  fields.  The  second  NOx 
model  (described  above)  was  used  as  the  postpro¬ 
cessor  to  the  pdf-predicted  temperature  and  chem¬ 
istry  fields.  Two  options  were  used  here,  one  with 


the  temperature  variance  calculations  and  the  sec¬ 
ond  one  with  the  mixture  fraction  variance  calcula¬ 
tions.  Table  2  reports  the  comparison  of  the  mea¬ 
sured  and  predicted  NOx  emissions.  While 
comparing  the  measurements  and  the  predictions  it 
should  be  remembered  that  a  large  part  of  the  fur¬ 
nace  exit  NOx  was  coming  from  the  comburent  that 
contained  110  ppmvd  NOx.  Relatively  little  net  NOx 
generation  (0.8  mg/s)  was  observed  in  the  measure¬ 
ments  if  compared  to  the  amount  (24  mg/s)  of  NO 
that  entered  the  furnace  with  the  comburent.  All  the 
NOx  submodels  performed  well,  and  the  predicted 
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Fig.  4.  Measured  and  predicted  NOx  concentrations  ppm,  vol,  dry.  The  predictions  marked  as  Fluent5-pdfll  case 
were  obtained  using  the  equation  for  temperature  variance  transport  (see  Table  2). 


net  NOx  generation  was  very  much  in  agreement 
with  the  measurements.  An  additional  computa¬ 
tional  run,  using  the  IFRF-EDC  model  in  combi¬ 
nation  with  the  NOx)  submodel  described  above, 
was  made  with  the  inlet  (comburent)  NOx  level  set 
to  nil.  The  run  resulted  in  28  ppm  NOx  at  the  fur¬ 
nace  exit. 

Figure  4  shows  the  measured  and  predicted  in¬ 
furnace  NOx  concentrations.  The  peak  in  the  NOx 
concentration  was  at  the  boundary  of  the  air  and 
natural  gas  jet  where  the  highest  temperature  were 
both  measured  and  predicted.  The  predictions  were 
able  to  capture  well  this  behavior.  The  NOx  concen¬ 
tration  showed  a  peak  between  the  air  jet  and  the 
fuel  jet  at  the  two  first  traverses.  The  source  of  the 
NOx  (net)  formation  rate  was  located  in  this  thin, 
elongated  region  (flamelet)  between  the  air  and  the 
fuel  jet.  The  NO  net  formation  rates  were  several 
orders  of  magnitude  lower  than  in  conventional 
burner  technology  [13].  Using  the  EBU  and  the 
EDC  combustion  models  the  influence  of  the  NOx- 
reburning  mechanism  was  investigated.  The  NOx  re¬ 
duction  rates  were  4-5  orders  of  magnitude  lower 
than  the  NOx  formation  rates.  Additional  computa¬ 
tional  runs  were  made  with  the  thermal-NO  and 
prompt-NO  submodels  deactivated  (see  Table  2) 
showing  that  around  95%  of  the  furnace  exit  NO  was 
formed  via  the  thermal  path. 

The  Fuel  Jet 

We  have  already  observed  substantial  differences 
between  the  predicted  and  the  measured  tempera¬ 
tures  (Fig.  3)  in  the  fuel  jet,  irrespectively  of  the 


combustion  model  used.  Furthermore  these  differ¬ 
ences  seemed  to  be  of  little  importance  as  far  as  the 
NOx  predictions  were  concerned.  The  applications 
of  more  sophisticated  turbulence  models  (see  above) 
brought  little  changes,  and  no  substantial  changes 
were  observed  when  the  kinetic  turbulence  energy 
at  the  inlet  of  the  fuel  jet  was  increased  by  a  factor 
of  three.  This  was  so,  since  the  inlet  kinetic  energy 
was  rapidly  dissipated  before  the  fuel  jet  entered  the 
furnace  (the  computations  originated  at  100  mm  up¬ 
stream  of  the  furnace  front  wall;  see  Fig.  1).  We 
should  recall  that  our  computations  were  carried  out 
under  steady-state  conditions,  thus,  as  if  the  (weak) 
fuel  jets  were  stiff  and  straight.  The  laser  sheet  vi¬ 
sualization  applied  during  the  measurements  show 
the  fuel  jets  dancing  around  while  the  central  strong 
(comburent)  jet  remained  straight  and  stiff.  This  ob¬ 
servation  would  suggest  that  unless  these  non-sta- 
tionary  effects  are  accounted  for,  there  would  be  a 
little  chance  of  success  in  predicting  the  fuel  jet 
properties.  Figs.  2  and  3  show  an  additional  com¬ 
putational  run  (marked  IFRF-EDC-more-turb)  that 
predicted  both  the  temperature  and  the  methane 
concentrations  better  than  in  the  previously  analyzed 
runs.  In  this  run,  we  added  an  extra  source  term  in 
the  transport  equation  for  the  kinetic  energy  of  tur¬ 
bulence.  The  source  was  active  only  in  the  fuel 
(methane)  jet.  This  is  indeed  an  ad  hoc  modification 
and  as  such  is  unjustified  in  the  steady-state  calcu¬ 
lations.  In  our  opinion,  there  are  strong  indications 
that  modeling  of  the  strong  jet-weak  jet  interactions 
would  require  a  model  that  should  account  for  the 
enhanced  entrainment  of  oxygen  and  enhanced  rate 
of  fuel  combustion  due  to  the  non-stationary  behav¬ 
ior  of  the  weak  (methane)  jet.  Interestingly  enough, 
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running  the  NOx  postprocessor,  after  improving  the 
temperature  and  methane  concentration  predictions 
in  the  weak  (methane)  jet,  did  not  result  in  changes 
in  the  predicted  furnace  exit  NO  emissions,  as  re¬ 
ported  in  Table  2.  The  NO  source,  in  terms  of  its 
location  and  magnitude,  remained  unaltered  since 
the  peak  in-furnace  temperatures  remained  the 
same. 


Conclusions 

Steady-state  numerical  simulations  of  a  0.58  MW 
burner  that  operated  in  a  flameless  oxidation  mode 
have  been  carried  out.  Three  combustion  models 
have  been  tested: 

•  Eddy-break-up  model  with  a  two-step  reaction 
scheme 

•  Eddy  dissipation  concept  model  with  chemical 
equilibrium 

•  pdf/mixture  fraction  model  with  equilibrium,  non- 
adiabatic  look-up  tables 

The  nitric  oxide  formation  has  been  investigated 
with  two  different  NOx  models  both  considering 
thermal  and  prompt  NOx  formation.  In  one  of  the 
NOx  models  the  NO-reburning  mechanism  has  been 
included.  The  model  results  have  been  compared 
both  with  the  in-furnace  measurements  (including 
temperature,  gaseous  species,  velocity,  radiative  heat 
flux)  and  with  the  furnace  exit  parameters.  The  fol¬ 
lowing  has  been  concluded: 

1.  With  the  exception  of  the  small  region  located 
within  the  natural  gas  jet,  the  computations  have 
resulted  in  predictions  of  good  quality.  The  com¬ 
puted  NOx  emissions  have  been  in  very  good 
agreement  with  the  measured  values  for  a  num¬ 
ber  of  combustion-model/NOx  postprocessor 
combinations. 

2.  The  NOx  is  formed  in  an  elongated,  thin  shear 
layer  located  between  the  natural  gas  jet  and  the 
air  jet.  To  obtain  good  quality  NOx  emissions  pre¬ 
dictions,  it  was  enough  to  predict  correctly  the 
time-mean  temperature  and  oxygen  concentra¬ 
tions  in  this  (flamelet)  region.  The  subsequent  ap¬ 
plication  of  the  f- pdf  function  either  to  the  tem¬ 
perature  or  to  the  mixture  fraction  resulted  in 
correctly  predicted  NOx  emissions. 

3.  The  NOx  is  formed  mainly  by  the  thermal  path. 
The  prompt  NO  is  of  little  importance.  The  NO- 
reburning  mechanism  is  negligible  in  the  ana¬ 
lyzed  case, 

4.  There  are  strong  indications  that  non-stationary 
behavior  of  the  weak  (fuel)  jet  must  be  accounted 
for  to  improve  the  predictions  within  the  jet. 
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COMMENTS 


Bassam  Dally,  University  of  Adelaide,  Australia.  I’d  like 
to  comment  that  I  measured  similar  flames  and  found  that 
finite-rate  chemistry  is  important  and  hence  the  model  you 
used  may  not  be  applicable.  It  seems  to  me  that  the  initial 
conditions  at  the  jet  exit  may  contribute  to  the  discrepancy. 
Can  you  comment  on  that?  Do  you  think  that  including 
radiation  in  your  model  will  improve  your  predictions? 


Authors  Reply.  Two  years  ago,  we  also  thought  that  the 
finite-rate  chemistry  might  be  important  [1].  As  for  today, 
after  carrying  out  further  computations  (Ref.  [12]  in  pa¬ 
per),  we  understand  that  the  finite-rate  chemistry  cannot 
be  responsible  for  such  large  differences  between  the  mea¬ 
sured  and  predicted  properties  of  the  natural  gas  jet.  The 


uncertainties  in  the  inlet  conditions  of  the  jet  are  not  re¬ 
sponsible  for  these  differences.  In  our  numerical  simula¬ 
tions,  we  use  the  ordinate  transfer  radiation  method.  The 
absorption  coefficients  are  calculated  spectrally  or  using  a 
gray  gas  approximation.  In  either  case,  the  radiation  alters 
the  temperatures  in  the  natural  gas  jet  by  30-50  °C,  while 
the  observed  discrepancies  are  substantially  larger.  As  we 
now  understand  the  problem,  the  non-stationary  behavior 
of  the  weak  natural  gas  jets  is  responsible  for  the  discrep¬ 
ancies. 
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